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Abstract. The hypothesis that AGN jets might be the sources of the ultra-high energy cosmic 
rays has originally motivated the venture of Te V gamma ray astronomy. Surprisingly, after the 
discovery of TeV emission from blazars the attention has shifted to more traditional explana- 
tions which do not involve energetic hadrons, and there is even common believe that a hadronic 
interpretation is disfavored by observations. It is shown here that this is not the case, and that 
the currently observed spectra and variability features of blazars can be perfectly understood 
within hadronic blazar models. I also discuss how hadronic models might be observationally 
J> ■ distinguished from common leptonic models, and point out some interesting aspects which 

04 ' could be relevant for the understanding of the differences between blazar classes. 
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O : WHY HADRONIC MODELS? 
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AGN jets and the origin of cosmic rays 

Cosmic rays are observed up to the enormous energy of 3 x 10 20 eV, and both theoret- 
ical and observational arguments suggest an extragalactic origin of these most energetic 
particles [1]. Since the cosmic ray arrival directions are largely randomized by Galactic 
and extragalactic magnetic fields, their sources cannot be easily identified from direct 
observations. However, if we assume that they gain their energy by acceleration (rather 
than by quantum processes, i.e., the decay of superheavy particles) their extreme maxi- 
mum energy allows to derive quite restrictive selection criteria for possible acceleration 
sites. Most acceleration scenarios in astrophysics, for example Fermi acceleration [2], 
assume that particles are magnetically confined for some time t acc in the accelerating 
region. This implies three fundamental constraints on the maximum proton energy, 

E cr < eBRT , (la) 
E CI < \m v c 2 Y{m v j m e ) ^ r]B c /a { B (lb) 
E CI < ijeBXczxT, (lc) 

which we call (a) the confinement limit, (b) the synchrotron limit, and (c) the GZK limit. 
The confinement limit simply expresses that the particle gyro-radius in the magnetic 
field B is smaller than the system size R. The synchrotron limit expresses that t acc is 
smaller than the synchrotron loss time, where B c « 4xl0 13 G is the critical magnetic 



field, and a f 1/137 the fine structure constant. The parameter 77 = AE/E is the 
fractional energy gain in a Larmor time, t L = E/eBc, and one can show that for Fermi 
acceleration generally 77 < 1 (see [3, App. D] for some discussion and references). The 
GZK limit expresses that proton acceleration must be faster than photohadronic losses at 
the microwave background, called the Greisen-Zatsepin-Kuzmin or GZK effect [4], and 
Agzk ~ 10 Mpc is the appropriate attenuation length of protons with E > 10 20 eV [5]. 
The GZK limit dominates over the confinement limit for sources on supercluster scales 
(>10 Mpc), but it applies also to somewhat smaller scales (e.g. clusters of galaxies) since 
in such objects usually rj <C 1 [6]. 

Only very few cosmic sources have been found which may satisfy all three conditions 
for E CT > 10 20 eV, and all of them are connected to strong shocks in relativistic flows 
(r ~ 10—1000): (a) the termination shocks of extended jets in radio galaxies [7], (b) 
compact jets in blazars [8,9], and (c) internal or external shocks in fireballs proposed 
to induce Gamma-Ray Bursts [10]. All of them are known as powerful emitters of non- 
thermal photons, in particular gamma rays. If we adopt the assumption that there is a uni- 
versal ratio between cosmic ray and non-thermal photon emission (which is rather sim- 
plistic, but not unreasonable as an estimate), we could use their total contribution to the 
observed extragalactic gamma rays as a scale of their total power as cosmic ray sources. 
This clearly favors AGN over GRB, because AGN are known to produce 10—30% of 
all extragalactic gamma-rays (counting both resolved sources and diffuse background), 
while resolved GRB contribute less than 1%. Due to the GZK effect, another important 
selection criterion for UHECR sources is proximity. This disfavors powerful quasars, 
and also GRB if their cosmological distribution follows the star forming activity, since 
in both cases most power is emitted at large redshifts. In contrast, BL Lac objects and 
their proposed unbeamed counterpart, FR-I radio galaxies [11], seem to be more fre- 
quent in the current epoch than at large redshifts [12]. Since their jets may also satisfy 
the condition for the acceleration of UHECR, and they contribute significantly to the 
non-thermal radiation in the universe, blazars and FR-I radio galaxies may be regarded 
the primary candidates for the origin of the highest energy cosmic rays. 

Hadronic vs. leptonic gamma-ray emission in blazars 

Energetic hadrons can lead to the emission of gamma rays via pp interactions with sur- 
rounding gas, or p-y interactions with ambient photons. This leads to the production of 
secondary e ± pairs, or mesons like 7r ± and 7r°, which eventually decay into e ± pairs, pho- 
tons and neutrinos. Electrons or pairs can produce high energy photons by synchrotron 
(or Compton) processes. The photons can either escape from the jet or produce new 
pairs in 77 — > e + e~ processes, which subsequently radiate a new generation of pho- 
tons. In particular synchrotron-pair cascades of this kind are important if interactions of 
UHECR protons in AGN jets are considered, where they shift the energy in secondary 
radiation down from the extreme proton energies to the observable gamma-ray regime. 
This mechanism has been coined PIC for "proton induced cascades" [13]. 

The term "hadronic blazar models" is used for a large variety of models for the gamma- 
ray production in blazars, not all of which involve UHE cosmic rays. All in common is 



just that they propose energetic protons as the main carrier of dissipated energy in the jet, 
rather than so-called "leptonic models", which assume the bulk of the energy available 
for radiation in electrons or pairs. Some hadronic models assume gamma-ray produc- 
tion in pp interactions invoked in the collision of jets with surrounding gas clouds, or in 
a very massive jet itself [14]. In a different class of models, protons have been suggested 
as being responsible just for the injection of energetic electrons, which then produce the 
observed photon emission by a synchrotron-self Compton (SSC) mechanism [15]. These 
models offer an interesting explanation for the observed strong variability in gamma-ray 
blazars due to intrinsic instabilities. They usually require only moderate proton energies, 
but very large densities of relativistic protons in the jet. 

Following the motivation given above, I want to focus here on a different kind of 
hadronic models in which UHE cosmic rays in the jet interact with low energy target 
photons. They are split into two classes: (a) target photons are produced by synchrotron- 
radiating electrons co-accelerated with the protons [8], and (b) external target photons 
are present in the vicinity of the jet, as for example thermal photons emitted from an 
accretion disk or a warm dust torus [16]. In a realistic scenario, photons from both 
sources would be present. The local ratio of the external and internal photon energy 
density in a reference frame comoving with the jet differs hereby from the the ratio of the 
observed external and internal luminosities by a factor r? et (the corresponding number 
ratio of external to internal photons scales with Tf et ). It is therefore hard to constrain 
the local comoving density of external photons in the jet by observations. However, 
an indirect constraint on the external photon density has been pointed out by Protheroe 
and Biermann [17]: using a common accretion disk/torus model to estimate the radiation 
fields in AGN, they show that the emission of TcV photons would be strongly suppressed 
by 77 absorption, if the emission region is close enough to the AGN that external photons 
could be relevant for pj interactions. This is inconsistent with the observed spectra in 
TeV blazars like Mrk501 or Mrk421 (see below). Nevertheless, external photons can 
still be important in high luminosity EGRET blazars (like 3C 279) [9]. 

For the discussion of TeV blazars, we can therefore concentrate on so-called 
synchrotron-self proton induced cascade (SS-PIC) models, in which UHE protons in- 
teract with synchrotron photons emitted by electrons accelerated in the same process as 
the protons. One clear consequence of such models is that the gamma-ray emission due 
to the PIC process is in competition with SSC emission of the electrons. The relation of 
the luminosities of both processes is expressed by 
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where u p = N p dN p /dE p is the energy density of protons at the maximum energy, E p , 
and M ph is the bolometric energy density of the target photons (i.e., the synchrotron pho- 
tons emitted by primary electrons), and B and R are magnetic field and size of the 
emission region, respectively. The relation gives in fact only an upper limit on the PIC 
contribution, since approximate equality assumes E p = eBR. We see that, if the UHE 
proton content of the jet is significant (u p (E p ) > it ph ), and the maximum energy suffi- 



cient to explain the observed UHECR {eBRY )Ct > 10 20 eV, assuming a jet Lorentz factor 
Fjet ~ 10), PIC emission will dominate over SSC at compact jet scales (R < 10 16 cm) 
— and vice versa. For the likely case of adiabatic scaling of a dominantly transversal 
magnetic field with the jet radius, B oc it! -1 , the relative contribution of SSC emission 
would increase on larger scales. Therefore both processes may contribute to the observed 
radiation, although emission which is variable on short time scales (requiring small R) 
would tend to be dominated by PIC. Defining clear signatures to distinguish PIC from 
SSC emission can therefore make the hypothesis that blazars and FR-I radio galaxies are 
the sources of UHE cosmic rays testable by GeV-TeV gamma-ray observations. 



SPECTRAL PROPERTIES OF SS-PIC MODELS 

Synchrotron-pair cascades in power law target spectra 

At frequencies below the X-ray band, blazar spectra are well represented by broken 
power laws, in some cases even by single power laws down to the sub-mm regime [18]. 
These photons play a triple role in SS-PIC models: (i) they are targets for the initial 
P7 interactions, (ii) they are targets for the propagation of synchrotron pair cascades, 
and (iii) they allow clues on the particle spectrum of the primary accelerator. Let us, 
for simplicity, assume that the soft target photon number spectrum is described by a 
single power law, dN ph /dE oc a t ~ 1 is called the energy index of the target 

spectrum. The stationary electron spectrum producing these photons by synchrotron 
radiation must then have the form dN e /dE oc Q e ,acct sy n E~° le ~ 1 with an energy index 
a e = 2a t , where Q acc oc E~ a!iCC ~ l is the number of electrons injected by the accelerator 
per unit time, and t syn oc E~ l is the synchrotron cooling time of the electrons. This 
means that the spectrum injected by the accelerator is a power law with energy index 
«acc — a e — 1. In contrast to electrons, cooling of protons is usually dominated by 
adiabatic [3] or advection losses [8]. Hence, the proton cooling time t p can be considered 
as energy independent, and the stationary proton spectrum is dN p /dE oc Q Piacc (E). If 
the accelerated protons have the same spectrum as the electrons, the proton power law 
energy index is a p oc a acc = 2a t — 1, but we discuss also different choices of a p . 

The opacity for the photohadronic production is r P7 = t p /t pl , where t P1 oc E at 
is the according loss time scale. The major channel for hadronic gamma production 
is p7 — > 7T° + . . . with the subsequent decay n° — >■ 77. The resulting photons are 
dominantly absorbed by the soft target photons through 77 pair production. The 77 
opacity has the same energy dependence as r P7 for constant t p , that is r 77 oc E~ at , 
hence the stationary energy distribution of the primary photons from pion decay is 
dN!®/dE = dN p /dE(T 7J /r P7 ) oc dN p /dE. Saturated production of pairs in 77 col- 
lisions (i.e. t 77 3> 1) then leads to [19] 

diVf 1 r° , , f ,,dN® _ a± _, J x a«-ap-2 for a > a t 
— - — oc — / dx r 77 (x ) - ' oc x ± = < _ 2 c ^ (3) 

dx x 2 hx 77V 1 dx' \ x 2 for a p < a t 

where x is the energy of photons or pairs in units of m e c 2 , and x = ■y p m 7V /2m e ~ 1007 P 
is the maximum x of the primary injected photons. The pairs produce a new generation 



of photons by synchrotron radiation, which are distributed in the stationary, saturated 
case as 



oc — — — — oc x 1 
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i.e., as a power law with energy index a,\ = \a± + a t . Q^ n (x) oc :r~ a± / 2 is the injected 
number of synchrotron photons per unit time at energy m e c 2 x by the first generation of 
pairs. The stationary photons can inject a second generation of pairs, and so on. This 
synchrotron-pair cascade continues, and with each step the characteristic photon energy 
is reduced by 
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The latter condition expresses the classical limit of synchrotron radiation and is mostly 
fulfilled in hadronic blazar models; in the non-classical case the cascade propagates ap- 



proximately as x r 



\x n -i. Thus, the energy is rapidly reduced in each step once 



x <C B c j B ~ 10 , and since r 77 oc x at the opacity for the production of subsequent 
cascade generations quickly decreases. The cascade becomes unsaturated at a photon 
energy ~m e c 2 x 77 , defined by r 77 (:r 77 ) = 1, and dies out rapidly for x < x 77 . The 
photon spectrum emerging from the emitter is [19,13] 
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where the sum extends over all cascade generations and n* is the last generation with a 
highest photon energy x n * > rc 77 (x n is determined by repeated application of Eq. (5) 
on x ). It can be shown that for typical blazar conditions n* = 3—4, and that Q^*J 
dominates the emission around x 77 [20]. 

The term in brackets in Eq. (6) has the meaning of a mean escape probability of the 
photons from the emission region, and is strictly correct only for a plane-parallel ge- 
ometry. However, the asymptotic behavior of the spectrum would be the same in any 
geometry, that is, x 77 marks a break in the spectrum by Aa = a t , which we call the 
opacity break. This broken power law shape is typical for photospheric emission, which 
is in the nature of the SS-PIC models where the soft photons are responsible both for 
the production and absorption of gamma rays. It is directly observable unless there is 
significant absorption by external photons surrounding the jet, which would lead to an 
exponential cutoff oc exp(— r 77jCxt ). The observed GeV— TeV spectra of TeV-blazars 
support in most cases a broken power law shape, as expected from the SS-PIC model 
for m e c 2 x^Tj Ct ~ 1 TeV, but are inconsistent with a rapid cutoff in this regime [18]. 
This has been used as an argument against external PIC models for TeV blazars [17]. An 
observed opacity break in the TeV regime is also expected from a direct determination 
of t 77 for typical blazar parameters [8], and from general considerations concerning the 
efficiency of hadronic blazar models [9]. 



Robust features ofhadronic blazar spectra 



If we ignore the details of the cascade spectra in the vicinity of cutoffs or breaks, the 
spectral propagation can be described in terms of a simple algebra. If we call a n (x n ) the 
local spectral index in the vicinity i ~ i n , we can introduce the relations 

) = /+[a n -i(^n-i)] for x n > x 77 (7a) 

) = f-[a n -i(x n -i)} for x n < x 77 (7b) 

with 

f±[a n ] = max{|(a n ±a t + 1),|} , (7c) 

where the dimensionless photon energy x n propagates following Eq. (5). The spectrum 
of the /c-th cascade is then given by /+[<%] above the opacity break at x 77 , and by 
/-[/+ _1 [ a p]] below (but above x'^ = Bx^ /4B C <C x 77 ), where /£ denotes the A; -fold 
iterative application of / + . 

If we ask which cascade generation is dominant at a given energy we encounter a 
problem related to the energy propagation equation (5): for any given primary spectral 
feature at x , the position of its "image" in the ra-th cascade is dependent on the magnetic 
field in the emitter by x n oc B 2 " -1 . This means, for example, that the cutoff of the 
4-th cascade generation is x A oc B 15 , so that a variation of B by a factor of 2 over the 
emitting region would "smear out" the value of x A by more than 4 orders of magnitude. 
Fortunately, it is still possible to make some quite robust predictions on hadronic blazar 
spectra, since the cascade generation spectra converge quickly with k to /£° [a p \ = l+a t , 
and /_ [/2°[a!p]] = 1, independent of a p . Moreover, the power contained in each cascade 
generation is approximately equal for n < n* + 1 and rapidly decreases for larger n, so 
that the spectrum below x 77 is not strongly changed by adding generations with n > n*. 
This allows the statement that hadronic gamma-ray blazar spectra are described by a 
broken power law with energy indices a ~ 1 below, and a ~ 1 + at above a break 
observed at ~1 TeV, where a t is the IR-X energy index of the source. 

So far we did not consider the cutoff in the target photon spectrum, x t , which is gen- 
erally observed at an energy m e c 2 x t Tj et between in the optical and X-ray regimes. We 
consider the target photon spectrum for x t > x ( as a power law with a steep energy 
index, a' t > a t + 1. Gamma rays with x n < x^ 1 then induce cascade photons below 
x = B/(AB c xf) with an index a' n+1 R3 f'-[f+[<^p]] = \, where f'_ is the cascade opera- 
tor defined in Eq. (7c) replacing a t with a' t (for x < x 77 ). Since this result in independent 
of the detailed value of a' t as long as a' t > a t + 1, we can consider it as valid for any steep 
cutoff, e.g., an exponential one. The energy m e c 2 xTj et marks a break in the spectrum, 
which is observable above the cutoff of the primary electron emission (x > x t ) if 

e t = m e c 2 T sct x t < 100 eV 

For B < 100 G this transition corresponds to the observationally motivated distinc- 
tion between so-called low-energy-cutoff blazars (LBLs) and high-energy-cutoff blazars 
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FIGURE 1. Generic SS-PIC spectra for HBL (left) and LBL (right) type blazars. Dashed high-energy 
lines show a typical spectral shape of an SSC model for comparison. 



(HBLs) [11]. In the SS-PIC models we therefore expect fundamentally different shapes 
for the high energy spectra of these classes, as illustrated in Figure 1 . Obviously, the 
spectrum for LBLs can easily be confused with an SSC spectrum, while a clear differ- 
ence is predicted for HBLs in the MeV gamma-ray regime. 



Narrow cascades: proton and muon synchrotron radiation 

Apart from n° decay, energetic photons can also be produced by synchrotron radiation 
of the protons. Moreover, it has been shown that also the most energetic muons produced 
by the decay of photohadronically produced charged mesons can lose a significant frac- 
tion of their energy in synchrotron radiation in magnetic fields typical for SS-PIC blazar 
models, prior to their decay [3]. The total power ratio of synchrotron radiation compared 
to the 7T° cascade injected by protons of the energy Ecr in the observer's frame is 
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where L IR is the observed infrared luminosity of the blazar. For muon synchrotron radi- 
ation the ratio is L^ sya /L^o « 2 because of the photohadronic branching ratios [21]. 

If the jet contains cosmic ray protons up to an energy E CT ~ 3xl0 20 eV in the ob- 
server's frame, B > 10 G and the typical parameters assumed in Eq. (9) otherwise, pro- 
ton synchrotron radiation extends up to observable energies of Bx 2 p m z J?- jCt / (AB c m z p ) > 
300 GeV. Muon synchrotron radiation requires a minimal dimensionless energy x* = 
E^/m e c 2 Tj Ct ~ 5xlO n for B ~ 10 G [3], and therefore extends from 30 GeV to 3 TeV 
using the same parameters. Since we have argued above that the opacity break is in the 
same regime, a significant fraction of this radiation can be reprocessed in a further cas- 
cade generation and would appear in the hard X-ray regime up to energies of ~10 keV 
and ~100 keV for reprocessed proton and muon synchrotron radiation, respectively. 

The interesting aspect about these cascades becomes obvious when we consider their 
spectrum below the opacity break. Here, we assume for simplicity that proton syn- 
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FIGURE 2. Schematic spectra of narrow cascades superposing the broad band emission of ordinary ir° 
cascades and primary electron synchrotron radiation — Spectral indices correspond to ^i-induced cascades 
for a p — 1, arrows indicate the cascading of spectral features. The proton-synchrotron cascade is omitted 
for clarity, but has a similar structure at photon energies about a factor of 10 lower (see also [23]). 



chrotron radiation is significant, but not dominant over adiabatic or advection losses. In 
this case, the proton synchrotron spectrum has an energy index a Pt i = \a p ~ 0.5. For 
muons synchrotron cooling is always dominant (otherwise it is suppressed due to muon 
decay), and we obtain = f~[a p ] ~ 0.5 for x 77 > x > x\ = Bx* 2 ml/ '(4£? c m^), and 
a^i = — | below from synchrotron emission below the characteristic frequency. Hence, 
both processes produce spectra very much flatter than those arising from ordinary, it 
induced cascades. Above the opacity break, the spectra steepen by a t ~ 1, which means 
that the peak of the luminosity is reached at the opacity break. Here, the dominance 
of this component over the broad it cascades is even stronger than expected from the 
power ratio factors derived above, since it is first generation synchrotron emission, while 
the 7T° emission from the same protons is broadened through the cascading process. The 
spectral indices expected from the narrow cascades below and above the opacity break 
(Fig. 2) are consistent with the typical EGRET and TeV indices, respectively, observed in 
Mrk421 and Mrk501 [18]. The spectral indices of the second cascade are obtained from 
a second application of the operator /~, yielding again values around 0.5 steepening by 
\a t above an X-ray break observable at m e c 2 Tj Ct x'^ with x'^ « Bx 2 ^/AB C (see Fig. 2). 
These compare well to the indices observed by Beppo-SAX in three prominent flares of 
Mrk501 in April 1997 [22]. The luminosity in the X-ray peak depends sensitively on 
the energy of the opacity break, which determines how much energy of the TeV peak 
is reprocessed. For example, in the SS-PIC model by Miicke and Protheroe [23] the 
opacity break is at m e c 2 x^Tj et ~ 25 TeV, which explains why their X-ray peak from 



reprocessed proton synchrotron radiation is strongly suppressed. For suitable parame- 
ters, however, narrow PIC emission can produce a two-bump spectrum with comparable 
peaks the TeV and the X-ray regime, as illustrated in Figure 2. 

VARIABILITY AND CORRELATED FLARES 

Flares and the quiescent background 

A major result of the simultaneous multiwaveband observations of Mrk421 and 
Mrk501 was that their variability in the X-ray and TeV-regime is largely correlated, and 
stronger than in other wavebands. This has been used as an argument against a hadronic 
interpretation of their gamma-ray emission, since ordinary cascade models expect a quite 
model independent spectral index (as explained above), which implies that the gamma- 
ray variability of such sources should be comparable at all gamma-ray energies. 

Obviously, this picture changes if we consider narrow cascades. In order to play a 
dominant role in the emission, these require conditions which allow to produce cosmic 
rays of extremely high energies. Let us assume now that such conditions are not always 
present in the jet, but only in some confined regions for a limited time. Of course, 
we may still assume that also in other regions the jet accelerates protons and electrons, 
albeit not to such high energies. This would lead to a permanent "glow" of the jet, 
which is dominated by the emission of primary electrons and n° induced cascades with a 
spectrum comparable to the case illustrated in Fig. 1 for HBLs, which is consistent with 
observations of the quiescent emission of Mrk421. A short flare which contains UHE 
protons able to produce narrow cascades, and which is energetic enough to compete with 
the total emission of the rest of the jet, would then cause the strongest variability in the 
energy regimes where the flare spectrum peaks, and these are the X-ray and TeV regimes. 
At other wavebands, the variability may be low, or not present at all since the emission 
there may continue to be dominated by the quiescent background (see Fig. 3). 

Opacity and variability 

The assumption of the existence of an opacity break in or below the TeV regime is 
vital for a successful explanation of blazar spectra in an SS-PIC model. This has some 
very interesting implications if the emission is variable. Let us assume a nearly plane- 
parallel geometry of length R for a region emitting photons simultaneously over its entire 
volume, over a time scale t rad <C Rj c. In the optically thin case, this would induce a flare 
of duration R/c owing to the run-time differences of photons. If the emitter is opaque 
at some some energy x, i.e., t J7 (x) 3> I, only photons emitted within a distance i?/r 77 
can reach the observer, causing a flare of duration i?/cr 77 . 

Applying this to SS-PIC models with an opacity break below TeV, we would there- 
fore expect that the observed multi-TeV variability becomes systematically faster with 
increasing photon energy, while the variability time scales are energy independent at 
X-ray or EGRET energies that are below the opacity break. Of course, the simple rela- 
tion T var = T var /r 77 for r 77 > 1, and T var = T var for r 77 < 1 is unlikely to apply 




photon energy 

FIGURE 3. Correlated variability on blazars as a result of flaring, narrow PIC emission superposing a 
non-peaked background from the surrounding jet. 

to SS-PIC, since both cascade injection and propagation involves time scales ~i?/c. 
Nonetheless, the effect remains qualitatively in the sense that it reduces the variability 
time from ~ t rad + R/c to ~ t rad for r 77 3> 1. This may explain some recent simul- 
taneous X-ray /TeV observations of Mrk501 which indicate a TeV variability about a 
factor of 3 shorter than at three different X-ray energies, which all show the same time 
curve [18]. In contrast, optically thin SSC models would expect that there is some X-ray 
energy which shows the same variability as observed in the TeV. 

SUPPLEMENTARY ASPECTS 

Particle acceleration in the SS-PIC model 

Some of the properties of the SS-PIC model discussed above on a purely heuristical 
base (i.e., the assumption that UHECR are produced) can be put on a more rigid foun- 
dation if we include the acceleration process. For jets the most reasonable assumption 
is Fermi acceleration at shocks or plasma turbulence [2]. If a particle with mass m is 
accelerated on a time scale t acc = rj'Hi,, where r\ has the same meaning as in Eqs. (1), 
then the cutoff of the synchrotron radiation emitted by the particle satisfies the condition 

e < nr] mc 2 /ctf ~ 400?7 mc 2 , (10) 

which is reached if the particle energy limited by synchrotron cooling (t acc = t syn ). 
Applying Eq. (10) to synchrotron radiating protons we find a maximum observed pho- 
ton energy of « [4 TeV]rj p (assuming r jet ~ 10). For Fermi acceleration we can write 
r/(E) ~ [P 5B(r L )/B] 2 , where (3 is the velocity of the shock or plasma wave, and 
5B(r L )/B is the fractional magnetic field turbulence on the scale of the particle gyro- 
radius, r L = E/eB (see [3, App.D], and references therein). Obviously, in order 



to explain photons observed above 300 GeV by proton synchrotron radiation we need 
r](E CT ) > 0.1, which requires relativistic shocks and strong turbulence on the largest 
scales in the system. In this picture, correlated flares in blazars are due to the appearance 
of transient relativistic shocks in a turbulent flow, while the background emission is due 
to continuous acceleration at omnipresent weak shocks or plasma turbulence. 

Applying Eq. (10) to accelerated electrons, we find a maximum synchrotron photon 
energy of ~ [2GeV]?7 e for rj Ct ~ 10, which implies r] e ~ 10~ 9 — 10~ 5 to explain the 
observed IR-X cutoffs in blazars. The large difference between r\ v and r\ e can be un- 
derstood from the theory of plasma turbulence, noting that the electrons probe very 
much smaller scales of the turbulence than the protons. If the turbulence is described 
by r)(E) oc [5B(r L )/B] 2 oc E», we find 



e e,syn ~ e p,syn 



i (32/+2)/(2-y) 

m e /m p , (11) 



which is obviously independent of Doppler boosting. Biermann and Strittmatter [24] 
pointed out that the near infrared cutoffs observed in many quasars (e.g., LBLs) can be 
understood if they accelerate protons to the GZK limit, and the plasma turbulence is 



described by a Kolmogorov spectrum (y = p. In terms of Eq. (11), we would obtain 
e e ,syn ~ 1 eV for y = | and e P)Syn =10 GeV, which is consistent with the observations 
of LBLs and would imply i] p ~ 10~ 3 . The Kolmogorov spectrum applies to fully devel- 
oped hydrodynamical turbulence if the magnetic field does not significantly contribute 
to the total energy density of the fluid. In a magnetically dominated plasma a Kraich- 
nan turbulence spectrum would be expected [24], that is y = |. Combining this with 
ep.syn ~ 300 GeV which we required for blazars with correlated X-ray/TeV variability, 
we find e e syn ~ 10 keV, which corresponds to X-ray cutoffs typically observed in HBLs 
(note that the flare emission up to >100keV in Mrk501 is explained by hadronic cas- 
cades in this model). This allows an interesting explanation of the physical difference 
between these blazar classes: while LBL jets are energetically dominated by the hydro- 
dynamic motion of the plasma, and involve only non-relativistic shocks and/or weak to 
moderate turbulence, HBL jets are strongly turbulent, magnetically dominated flows, in 
which also relativistic shocks occur. More aspects of particle acceleration in the SS-PIC 
model are discussed in [23]. 



Gamma rays, cosmic rays and neutrinos 

One prediction which is unique to hadronic blazar models is the production of en- 
ergetic neutrinos with a luminosity comparable to gamma-rays. The detection of such 
neutrinos, in particular if correlated with blazar flares, would thus be a "smoking gun" 
for the hadronic scenario; unfortunately, the neutrino fluxes from single blazer flares 
are that low that this can hardly be expected within the next decades [3]. However, we 
would expect to find diffuse VHE neutrino fluxes comparable to the the diffuse extra- 
galactic gamma ray background (DEGRB), if a considerable fraction of the extragalactic 
gamma rays are of hadronic origin. This would also imply that a significant fraction of 
the UHECR flux is produced by the same process [9]. To utilize these relations to decide 



the nature of gamma-ray emission in blazars, we need to determine which fraction to the 
DEGRB they contribute — which is again a task for gamma-ray astronomy. 

In conclusion, although VHE neutrino observations will be very important to clarify 
the origin of cosmic rays, I still see the currently better technical possibilities to decide 
this question on the side of gamma-ray astronomy. To do this by revealing the nature of 
gamma-ray emission from blazars will require (i) a complete coverage of the gamma-ray 
wavebands from MeV to multi-TeV energies with sufficient spectral resolution, together 
with more campaigns allowing truly simultaneous multifrequency monitoring, and (ii) 
a comprehensive discussion of the data in the view of possible leptonic and hadronic 
explanations. 
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